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Abstract 



We discuss the low-frequency absorption arising from quenched inhomo- 
geneity in the superfluid density p s of a model superconductor. Such inhomo- 
geneities may arise in a high-T c superconductor from a wide variety of sources, 
including quenched random disorder and static charge density waves such as 
stripes. Using standard classical methods for treating randomly inhomoge- 
neous media, we show that both mechanisms produce additional absorption 
at finite frequencies. For a two- fluid model with weak mean-square fluctu- 
ations ((5p s ) 2 ) in p s and a frequency-independent quasiparticle conductiv- 
ity, the extra absorption has oscillator strength proportional to the quantity 
((5p s ) 2 ) I p s , as observed in some experiments. Similar behavior is found in 
a two-fluid model with anticorrelated fluctuations in the superfluid and nor- 
mal fluid densities. The extra absorption typically occurs as a Lorentzian 
centered at zero frequency. We present simple model calculations for this ex- 
tra absorption under conditions of both weak and strong fluctuations. The 
relation between our results and other model calculations is briefly discussed. 

PACS numbers: 



Typeset using REVTgX 



I. INTRODUCTION 



Some high-T c superconductors, such as Bi 2 Sr2Ca 2 Cu 2 08+ ( 5, absorb very strongly in the 
microwave regime even at temperatures T far below T c At a fixed frequency as a 

function of T, the real conductivity typically has two features: a sharp peak near T c , and 
a broad, frequency- dependent background extending to quite low T. The peak near T c is 
usually ascribed to critical fluctuations arising from the superconducting transition. 

The origin of the broad background is less clear. In conventional s-wave superconductors, 
there is no such background, because the material cannot absorb radiation at frequencies 
below the energy gap for pair excitations, 2 A. But in high-T c materials, which are thought 
to have a d x 2_ y 2 order parameter the gap vanishes in certain nodal directions in k- 

space. Hence, gapless nodal quasiparticles can be excited at arbitrarily low T, and these 
can absorb at low frequencies. But the observed absorption appears to be stronger than 
expected from the quasiparticles alone in a two-fluid model [Q. 

Several authors have suggested quenched inhomogeneities as the origin of this extra ab- 
sorption. Such inhomogeneities could arise from statistical fluctuations in the local densities 
of holes or impurities, or from the presence of superfluid-suppressing impurities such as Zn. 
Corson et al |j] found that such fluctuations in Bi 2 Sr 2 Ca 2 Cu 2 08+<5displace about 30% of the 
spectral weight from the condensate to the conductivity at finite frequencies. In addition, 
several explicit models have been presented which produce such a displacement. For exam- 
ple, Van der Marel and Tsvetkov || have calculated the fluctuation conductivity due to a 
periodic one-dimensional modulation of the phase stiffness. The present authors have 
calculated the fluctuation conductivity of a two-dimensional array of Josephson junctions 
with quenched randomness in the critical currents, using the Kubo formalism ||. They 
found an excess contribution to the conductivity below T c arising from this inhomogeneity, 
which arose from displacement of some of the superfluid contribution from zero to finite 
frequencies. 

More recently, Orenstein |J considered a two-fluid model of superconductivity with 
quenched, anticorrelated inhomogeneity in both the normal and superfluid densities; he 
found that some of the spectral weight of the conductivity is displaced from zero to finite 



frequencies, and gave explicit expressions for this extra conductivity. Han []T0]] has given a 
field-theoretic treatment of a randomly inhomogeneous two-dimensional d-wave supercon- 
ductor, using a replica formalism to treat the assumed weak disorder. He obtained an extra 
Lorentzian peak in the real conductivity, centered at zero frequency, arising from this dis- 
order, and a corresponding reduction in the superfluid density. All such models are made 
more plausible by imaging experiments which have directly observed spatial fluctuations in 
the local superconducting energy gap in underdoped cuprate superconductors |II . 



In this paper, we study the low-frequency complex conductivity of several models for 
two-dimensional superconductors with quenched inhomogeneities. We use a straightforward 
approach to obtain the effective complex conductivity, based on standard treatments of 
randomly inhomogeneous media in the quasistatic regime [E3]. From this approach, we 



obtain sum rules for the complex conductivity. We also find, in a number of cases, explicit 
expressions for the extra conductivity Sa e (u) due to the quenched inhomogeneities. Where 
the models are comparable, our explicit form for 5a e (uj) reduces to that of Han, even though 
obtained using a quite different approach. 
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The remainder of this paper is organized as follows. In Section II, we develop sum 
rules for the optical conductivity of a two-dimensional superconductor with quenched in- 
homogeneities, and apply them to several explicit models. In Section III, we present some 
approximate calculations for the effective conductivity of such superconductors |L3| . A con- 
cluding discussion is given in Section IV. 



II. CONDUCTIVITY SUM RULES 

A. Kramers-Kronig Relations 

We consider a two-dimensional superconductor, in which the local conductivity has nor- 
mal and superfluid contributions which act in parallel, but, in contrast to the usual two-fluid 
model, these terms are spatially varying. The complex local conductivity <r(x, uj) is taken 
as a scalar function of position x and frequency u) of the following form: 

a{u>,x) = a n {u,x) + (1) 

UJ 

a n might be a contribution from the nodal quasiparticles expected in a superconductor with 
d x 2_ y 2 order parameter symmetry, while p s represents the perfect-conductivity response of 
the superconductor. Our goal is to calculate the position-independent effective complex 
conductivity of this medium, which we denote 

a e (u) = a e>1 (uj) + ia e>2 (uj). (2) 

We first obtain a sum rule for a e> i. We start from the Kramers-Kronig relations obeyed 
by any function, such as a e (uj), which describes a causal response: 

cjV e , 2 (u/) 



<7 e ,i(u;) = -P / — -duj + doc (3) 

7T JO L) U — U) 1 



and 



a e 2 = en - ^ P r ff *w-*~ du/. (4) 

uj n Jo UJ' 1 — uj 2 

Here P denotes the principal part, and = Lim a; ^ 00 cr ei i(co'). We have assumed that 
Lim w ^ 00 cr n (co') position-independent, and have used the fact that, at high frequencies, a e ^ — > 
ip Sje /u, where p Sje is the effective superfluid density. At sufficiently large uj, the right hand 
side of eq. (|]) takes the form 

1 / 2 f°° \ 

0~e,2 -> - Ps,e + ~ We,^') - G^du' . (5) 

UJ \ 71 JO / 

We first write down a perturbation result for cr e (u;) , which will be needed below. If the 
spatial fluctuations in <r(u;,x) are small in magnitude compared to it spatial average, a av , 
then to second order in the |5cr|/|a ai ,| = \a(u,x) — a av (uj)\/\a av (uj)\, 

1 ((H 2 ) 



2 a n 



(6) 



where (...) denotes a space average of the quantity in brackets |JT2 
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B. Frequency-Independent Normal Conductivity 



Suppose first that a n (u) is real, nonzero, frequency-independent, and spatially uniform. 
Then = o~ n , and, at high u, however large the fluctuations in p s (x), the fluctuations 
in the local complex conductivity, a n + ip s (x)/u, are small compared to the space-averaged 
conductivity a av {uj) = a n + ip s<av /uj] we can thus use eq. (^). Since only p s is fluctuating, 
and not a n , eq. @ becomes 

a e K,a n ^ h - -= — , (7) 

to 2 \u 2 a n + iu>p S;av J 

where 5p s (x) = p s (x) — p av . Thus, to leading order in l/u, a e ^ ~ p s ,av/u. Equating this 
expression to the right-hand side of eq. fl5|), we obtain 

7T 

L cr el (u ) - <r n ] do; = - (p a ao - p SjE ) (8) 
o 2 

One consequence of eq. @ is that, if p s (x) is spatially varying, p s (x), p Sja „ > p s>e , the right- 
hand side of eq. (§) is positive, and i/iere to// 6e an additional contribution to a e>1 (u), above 

The sum rule (H) requires only that the conductivity fluctuations are asymptotically small 
at large u, and not necessarily that p s have small fluctuations. If, however, the fluctuations 
in p s are, in fact, small, then p s e can be calculated approximately using the analog of eq. 
(|), namely, 

Pv „ ^ _ I (JMT1 . (9) 

Ps,av 

From this approximate result, eq. fl8|) becomes 

*<(^) 2 > 



/■oo . 7T / 1 



/> 2 

rs,av 



(10) 



For fixed mean-square fluctuations ((Sp s ) 2 ) / ' p 2 s av , this integral is proportional to the average 
superfluid density p s ,av That is, the extra integrated fluctuation contribution to cr e i, above 
the mean contribution of the normal fluid, is proportional to p S)Cro . A similar result has been 
reported in experiments carried out over a range of mean superfluid densities PJ. 



C. Frequency-Dependent, Spatially Fluctuating Normal Conductivity 

Next, we consider a frequency-dependent cr n (u;,x), assuming that Lirn a ,_^. 00 cr n (a;, x) = 0, 
but allowing for a spatially varying o~ n (uJ, x). Then eq. (]4|) becomes 

p s , e 2to r°° <r e ,i(a/) A , 1 / 2 [°° \ 

CT ei2 (^) = P / ^ -> - p s , e + - / (T e ,l( u ) duJ ' I 11 ) 

uj -k Jo to' 2 — to 2 u \ n Jo J 

where the last result is again valid at sufficiently large to. We now assume a Drude form for 
cr„(u;,x), i. e., 
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p n (x)r n (x) 

<r n (a;,x) = - : -—, (12) 

1 - lUJT n (-X) 

where p n (x) is a suitable normal fluid density, and r n (x) is a quasiparticle relaxation time. 
Then at sufficiently high u>, 

i[ Ps (x) + p n (x)] 
<7 2 (o;,x) , (13) 

while o"i(u;, x) falls off at least as fast as 1/uo 2 . Similarly, a e ^ also falls off at least as fast as 
l/u> 2 at large uj. Hence, a e ~ io e ^ at high frequencies, where 

~ (, ,\ (P s P n ^ e fl A\ 

Ve,2[U) = • (14) 

UJ 

Here (p s + p n ) e denotes the effective superfluid density of a hypothetical material whose 
spatially varying local superfluid density is p s (x) + p n (x). Equating the coefficients of uj^ 1 
on right hand sides of eqs. ( |TT]) and (|T^) yields 

f°° IT 

/ a e ,l(uj)dUJ = - [(p s + p n )e - p Sj e] , (15) 

Jo 2 

which gives the integrated spectral weight of a e ^{uj). 

To calculate this weight, we first introduce <J n ,av = (cn)> an d note that 

Vn,l,av{u )du = ^ , (16) 

even if r n (x) is position-dependent. Hence, 

7T 

[(T e ,l(v') ~ <rn,l,av{u')\ duj' = - [(p s + p n ) e - p s , e - p n>a v] • (17) 

Thus <T e< i(u') again has some extra spectral weight beyond that of cr nj0 „(o;). 

We first estimate expression ( |TTD in the small- fluctuation regime where \Sp s \ <C p s ,av, 
assuming that p n is non-fluctuating. Following steps analogous to those leading up to eq. 
(|TCi|), we obtain 

OO 7T / 1 1 \ 

Ki(wO - cT nA , av (u')} dJ -> -((5 Ps ) 2 ) — . (18) 

4 \Ps,av Ps,av+PnJ 

In the limit p n ^> p S) av, this expression reduces to the right-hand side of eq. fllTD- In this 
regime, the extra spectral weight is again proportional to p s , a v A frequency-independent a n 
is a special case of this limiting behavior. In the opposite limit (p sa v ^ Pn), the behavior is 



[cTe,l,a.v(u') - Cr n>1 {u)')) dJ -> jp n ^ . (19) 

^ Ps 



Thus, in this limit, for fixed mean-square fluctuations in p s , the extra spectral weight is not 
proportional to p sm 
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If both p s and p n are spatially varying, eq. ( |17D can still be easily treated when spatial 
fluctuations in p s and p n are small. Then, to second order in these fluctuations, 

(Ps + Pn)e « Ps,av + Pn,a V ~ IMEI^D. , ( 20 ) 

^ Ps,av ~r Pn,av 

while p Sj6 is given approximately by eq. @. Substituting eqs. ( p0|) and (||) into eq. (|17D gives 



OO y]- 

We,i{u) - a n>1>av (uj))(ko = - 



i(8 Ps y) {[8(p s + Pn )fy _ (2l) 

Ps,av Ps,av Pn,av 



To interpret eq. fl2l|) , we assume that p s and p n are correlated according to p n = Xp s . 
Thus A = 1 and A = — 1 represent perfect correlation and perfect anticorrelation. The 
predictions of eq. (|2l|) are particularly striking in the latter postulated in some 

models of inhomogeneity @,[nj • In this case the sum p s + p n is not fluctuating. Then the 
second term on the right-hand side of eq. fl2~T|) vanishes, and the extra spectral weight due 
to the inhomogeneity is again proportional to {(Sp s ) 2 ) / p s ,av In the opposite case A = 1, the 
right-hand side of eq. may be negative, i. e. the total spectral weight is smaller than 
that given by o- n ,i,av 



D. Tensor Conductivity 



The preceding arguments apply equally well to a tensor conductivity, as expected if the 
superconducting layer contains charge stripes or other types of static charge density waves. 
In this case, the local conductivity should be a 2 x 2 tensor []T3] of the form a a/3 (u;,x) = 
J R- 1 (x)(T d (w) J R(x), where a d (u) is a diagonal 2x2 matrix with diagonal components a^i^ 1 ), 
ctb(o;), and -R(x) is a 2 x 2 rotation matrix describing the local orientation of the stripes. If 
this phase has a domain structure, the stripes would have either of two orientations relative 
to the layer crystalline axes, with equal probability. In a macroscopically isotropic layer, 
this domain structure leads to an scalar effective conductivity a e (uj). If, for example, the 
conductivity along the ] th principal axis is o~j(u) = cr n j(u) + ip s j/u, then the superfluid 
inhomogeneity produces an extra spectral weight at finite uj given by eq. ([17]). 



III. MODEL CALCULATIONS 

A. Weak Inhomogeneity 

We now supplement these sum rules with some explicit expressions for cr e (u;), beginning 
with weak inhomogeneity, |5cr| -C o av . In this case, a e is given by eq. (H). If only p s is 
spatially fluctuating, this equation reduces to eq. (0) and 

^)^a n + \{{5 Ps f) (22) 

1 Ps,av + U °n 

Thus, the extra absorption is a Lorentzian centered at u — 0, of half-width oc p s /cr n , and 
strength oc ((5p s ) 2 ). A result equivalent to eq. ( |22"|) has also been obtained by Han using 
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a replica formalism. Indeed, even the parameters of the two results are identical, provided 
we interpret ((5p s ) 2 )/2 in our model as equal to the quantity gA 2 /n in the model of Ref. 

Next, we consider the case where both p s and a n have weak spatial variations. We assume 
that o n is given by eq. (12), that p n is spatially varying, but not r n , and that p n and p s are 
correlated according to the rule 5p n = X5p s . Then a e is given approximately by eq. (6), with 

Pn,avTn , ^Ps,av /<->«-> \ 

o- a v = ~ : 1 (23) 

1 — lUJT n LO 

and 

((H 2 ) = ((^) 2 )f i - + T ^-) 2 - (24) 

Some representative plots of o e) \{uS) resulting from eqs. (H), (|23|) and ([24] ) are shown 
in Fig. 1 for A = 1 (perfect correlation) A = — 1 (perfect anticorrelation), and A = (no 
fluctuations in 8p n ). The case A = — 1 appears most plausible physically, since it implies 
that the sum p n + p s is non-fluctuating. It also produces the largest increase in cr e i at any 
given to. 

For ujr n <C 1, eqs. @, (|23|) , and fl24"D lead to the analytical result 

_ / \ ((^Ps) 2 ) cr n(W - 2Ar n p sa „ 

cr e i(u;J w er na „ H ■ — . (25) 

2 Ps,av + ^ 2 °lav 

Here a n ^ av = p n ,av T n and we have neglected corrections of order (ujr n ) 2 in the numerator. 

Eq. (|25|) is again identical with Han's result [ 1Q| , obtained using a replica formalism, if we 
make the identification ((5p s ) 2 )/2 gA 2 /7r. 



B. Temperature-Dependent Superfluid Density; Effective-Medium Approximation 

Next, we calculate a e {oj) using a simple model for the conductivity of a CuC>2 layer 
proposed in Ref. et al 0, but generalized to include inhomogeneities. We assume that 
cr(cj,x) is given by eq. (|l|), but for simplicity that it can have only two values, a a and Ob, 
with probabilities p and 1 — p. The corresponding conductivities are <Jj{oj) (j = A, B), where 

<Tj(u) = a n (uj) + i^ (26) 

Following Ref. |5[], we take cr n (u) = p n r n /(l — iuJT n ). We assume p n = aT, for T < T c , as 
expected for nodal quasiparticles in a d x 2_ y 2 superconductor, and p n (T > T c ) = aT c . We 
also take r" 1 = j3T. This linear T-dependence may be expected if l/r„ is determined pri- 
marily by quasiparticle-quasiparticle scattering in two dimensions ||16|| , rather than impurity 
scattering, which might give a T- independent l/r„. For the superfluid density, we write 
Ps,a = 7Ps,o, Ps,b = 7"Vs,o, where p B , (T) = p Sj0 (0)y / l - 2aT/p sfi {0). This form ensures 
that p Sj o decreases linearly with T at low T, as observed experimentally |T7[, and that it 
vanishes at a critical temperature T c as \JT C — T . We use the experimental values of the 
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parameters T c , <7o = p n T n , p s ,o(0), « and (3 (as quoted in Ref. 0) for our numerical estimates. 
Finally, we have arbitrarily assumed that the inhomogeneity parameter 7 = 3. 

We compute a e (u) using the Bruggeman effective- medium approximation (EMA) ||18|| : 

ct a (uj) + a e (u) ' cr B {uj) + a e (u) 

The physically relevant solution is obtained by requiring that cr e (u) be continuous in u and 
that a e> x(u) > 0. We also arbitrarily assume that p = 1/2. Our form for p sA and p S: B 
guarantees that the EMA solution p s , e {p = 1/2) = ^/p s ,APs,B = Ps,o- 

Fig. 2 shows the resulting a eA {u,T) for several frequencies ranging from 0.2 to 0.8 THz, 
the range measured in Ref. ||. Also shown are J^ max a e ^(uj,T)du! for u min /(27r) = 0.2THz 
and u max /(27i) = 0.8THz. Finally, we plot a n7 i(uj, T) for these frequencies, as well as 
Iu max CTn^dio. Clearly, a e ^(u) is substantially increased beyond the quasiparticle contri- 
bution, because of spatial fluctuations in p s . 



C. Percolation Effects 

We next consider cr l e (u;) near a percolative superconductor to normal (S-N) transi 
tion. Such a transition might occur, for example, in a single Cu02 layer doped by a non 
superconducting element such as Zn, as has already been discussed by several authors |TS| 



We assume cr(x) to be given by eq. (P, with a position-independent Drude a n given by 
eq. (0), and a p s (x) which has one of two values, p s or zero, with probabilities p and 
1 — p. We calculate a e {oj) using the EMA, eq. (p7|) . The EMA predicts an S-N transition 
at p = 0.5 = p c . For all p on either side of p c , u e ,i(p, oS) has a peak at uj = 0, whose height 
diverges and whose half- width decreases to zero, as p — > p c . This behavior is illustrated 
in Fig. 3, where we plot a Gi i(p,u) for several p > p c , as calculated in the EMA. The other 
parameters are indicated in the Figure caption. The behavior of a e> i is qualitatively similar 
to that shown in Fig. 2, but there is additional critical behavior near the S-N transition, 
which is absent from the weakly disordered system. 

More exact results near p c can be obtained with the help of a standard scaling hypothesis 
T2| , |2"0"|] . For the present model, it takes the form 



£ = '^(gN- (28) 

Here <7< and cr> denote the complex conductivities in the S and N regions of the layer; Ap = 
p — Pc s and t are standard percolation critical exponents, which depend on dimensionality, 
and possibly on other structural details of the layer; and F±(x) are characteristic scaling 



functions which apply above and below p c . The expected behavior [12] is F+(x) ~ 1 for 



1 < 1, and F-.(x) ~ x for x <C 1, while F±(x) ~ x <// ^ + *- ) for x 1. From these forms of the 



scaling functions, we can infer [^OJ that cr e ,i has a peak at to = whose half-width on both 



sides of the percolation threshold is 



|Aw| w —\Ap\ s+t , (29) 

0"n 
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which vanishes at the percolation threshold, where \Ap\ — > 0. The height of this zero- 
frequency peak diverges according to the law cr el (p, u = 0) — > \Ap\~ s on either side of the 
percolation threshold. For two-dimensional bond percolation, s = t m 1.30, while in the 
EMA in any dimension, s = t = 1. Thus, we expect that a ej i(p,u) is characterized near p c 
by a line of diverging height, and half- width which goes to zero near p c . Precisely this type 
of behavior is seen in the EMA curves of Fig. 3. 



D. Tensor Conductivity 



The results shown in Fig. 3 correspond to an in-plane scalar superffuid density p s (x) with 
a random spatial variation. If instead p s (x) were a 2 x 2 tensor with principal axes varying 
randomly with position (as might occur for random quenched domains of stripes), a similarly 
enhanced a e \ would be expected. If the two principal components of the conductivity 
tensor are given by eq. fl2"E|), and the stripes can point with equal probability along either 
crystallographic direction, then a e (u) is given in the EMA by eq. fl27|) with p = 0.5. The 
resulting a e) i{uj) would behave exactly like that shown in Fig. 2. 



IV. DISCUSSION 



It is useful to compare our results to those of previous workers, considering first the 
work of Han [ItJ. For weak fluctuations in p s , and a frequency-independent non-fluctuating 



a n , our results are identical to his, provided we identify our parameter {(Sp s ) 2 )/2 with his 
quantity gA 2 /ii. However, our formal derivation is more elementary than that of Ref. [ITU 



Because of this simplicity, we can easily consider a wide range of quenched inhomogeneities. 
Since the two approaches yield identical results for comparable types of inhomogeneities, 
we infer that our calculations (which are based on standard methods of treating random 
inhomogeneities in classical field equations such as those of electrostatics) can also be done 
using field-theoretic techniques. This field-theoretic approach may therefore be useful in 
treating other classical problems in heterogeneous media. 

We briefly discuss relation between the present model and that of Ref. J7J. In that 
earlier paper, the CuC>2 layer was treated as an array of small "grains" coupled together by 
overdamped resistively shunted Josephson junctions with Langevin thermal noise; a e (u,T) 
was then computed using the classical Kubo formalism. It consists of two parts: a sharp 
peak near the Kosterlitz-Thouless transition at T c , due to breaking of vortex-antivortex 
pairs, and a broader contribution, for T < T c , which exists only if the critical currents have 
quenched randomness. This second contribution corresponds to that considered here for 
spatially varying p s . For T well below T c , there are few vortices in the Josephson array; 
hence, the Josephson links in the model of Ref. |7| behave like inductors, and the array 
acts like an inhomogeneous LR network. This network is simply a discretized version of 
the two-fluid model considered here: the resistances R represent the normal quasiparticle 
channel, and the inductances L corresponds to the superfluid. 

To summarize, we have shown that, in a superconducting layer with a spatially varying 
p s , there is extra absorption beyond that from low- lying quasiparticles alone. Such quasipar- 
ticles are expected in a d x 2_ y 2 superconductor. For several models of weak inhomogeneity, 
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we have shown that the extra spectral weight at finite frequencies is proportional to the 
quantity ((5p s ) 2 ) / p s , as reported for T < T c in underdoped samples of BiSr 2 Ca2Cu20 8 + :r 
0. Furthermore, in this regime, the extra weight generally appears as a Lorentzian centered 
at zero frequency, with half-width proportional to p s /u n . Similar behavior is predicted for 
tensor inhomogeneities, as is expected in stripe geometries. We also predict that, near an 
inhomogeneous superconductor-normal transition, the height of the inhomogeneity peak in 
01 (u>) diverges while its half-width goes to zero. Thus, superfluid inhomogeneity may be the 
principal source of the extra absorption beyond the two-fluid model which is reported in the 
high-T c cuprates. 
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FIGURES 
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FIG. 1. Excess real part of effective conductivity, a e: \{u) — <Jn,i,av{u), when there are weak 
correlated spatial fluctuations in both p s and p n , normalized by the relative magnitude of such 
fluctuations. a n ,i,av = Pn,avT n /[l + ^ 2 Tn] is the normal-state contribution. We assume 8p n = X5p s , 
and show results for A = —1 (negative correlation, solid line), A = 1 (positive correlation, dotted 
line), and A = (no fluctuations in the normal fluid density, dashed line). The average superfluid 
and normal fluid densities, p s>av and p n ,av, are assumed equal; the conductivities are in units of the 
quasiparticle p n T n . 
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FIG. 2. (a) a eA (u;,T), for u/(2n) = 0.2 THz (solid line), 0.4 THz (dotted line), and 
0.8 THz (dashed line), for the model inhomogeneous superconductor described in the text. 
Also plotted are (b) E e ,i = ^'u^Hdw, (c) <7„,i(w,T), (d) E n ,i = <T n ,i(w)dw, (e) 

0- e ,i(u;,T)-c7 ni i(u;,T), and (f) S ej i-E nj i = f£™™[a e>1 (uj)-a ntl (uj)]du;, where w m i n /(27r) = 0.2Thz 
and L0 max /(2ir) = 0.8THz. 
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FIG. 3. Real part (T e i(p,u), in units of the quasiparticle p n T n , of the effective conductiv- 
ity of a two-dimensional composite of normal metal and superconductor, as calculated using the 
two-dimensional effective-medium approximation for several values of p > p c : p = 0.9 (solid line), 
p = 0.7 (dotted line) and p = 0.5 (dashed line). 
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